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ABSTRACT 


An approach to compute three dimensional flows using two stream functions is pre- 
sented. The method generates a boundary fitted grid as part of its solution. Commonly 
used two steps for computing the flow fields: (l) boundary fitted grid generation, and 
(2) solution of Navier-Stokes equations on the generated grid, are combined into a single 
step in the present approach. The method presented can be used to compute directly 3-D 
viscous flows, or the potential flow approximation of this method can be used to generate 
grid for other algorithms to compute 3-D viscous flows. 

The independent variables used are Xi a spatial coordinate, and £ and rj, values of 
stream functions along two sets of suitably chosen intersecting stream surfaces. The de- 
pendent variables used are the streamwise velocity, and two functions that describe the 
stream surfaces. Since for a three dimensional flow there is no unique way to define two 
sets of intersecting stream surfaces to cover the given flow, in the present study three dif- 
ferent types of two sets of intersecting stream surfaces are considered. First is presented 
the metric of the (x, £,»7) curvilinear coordinate system associated with each type. Next 
equations for the steady state transport of mass, momentum, and energy are presented in 
terms of the metric of the (x, £, t?) coordinate system. Also included are the inviscid and 
the parabolized approximations to the general transport equations. 
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I. INTRODUCTION 


Since the introduction of stream function by Lagrange (l) for two dimensional plane 
flows and by Stokes (2) for axisymmetric flows, the use of a stream function to study two 
dimensional flows has been extensive. Computation of two-dimensional incompressible po- 
tential flows using the stream function as the dependent variable and the space coordinates 
as the independent variables is well known and can be found in almost any introductory 
fluid mechanics book. The stream function has also been used in the computation of vis- 
cous flows. In the recent past, the stream function, along with the vorticity, has been 
used extensively to compute two-dimensional incompressible viscous flows. Patankar and 
Spalding (3) have used the stream function to construct the cross-stream coordinate for 
the computation of two-dimensional compressible ‘parabolic’ (boundary layer) type flows. 
Kwon and Pletcher (4) have used the stream function and the axial velocity as the depen- 
dent variables to compute two-dimensional incompressible separated channel flow. These 
are a few examples of the use of the stream function for the computation of viscous flows. 
More recently, streamlines of the incompressible potential flow corresponding to a given 
geometry have been used to construct boundary-fitted grid systems for the computation 
of viscous flows. The streamlines needed for the grid generation have been calculated by 
various methods. For example, Ghia et al. (5) generated the grid by the use of conformal 
mapping; Meyder (6), and Ferrel and Adamczyk (7), by solving the potential equation. A 
survey of the use of streamlines to generate a grid is included in the review article on grid 
generation by Thompson et al. (8). 

The corresponding development for three dimensional flows, that is the use of two 
stream functions to study three dimensional flows, so far has been limited. Several au- 
thors in the past have introduced two stream function to describe three dimensional flows. 
Among the pioneering works are the works of Clebsch (9), Prasil (10), Maeder and Wood 
(ll), and Yih (12). For the two dimensional plane flows it follows from the continuity 


equation 

du 1 du 2 
dx 1 dx 2 


( 1 . 1 ) 


that u l dx 2 — u 2 dx l is an exact differential of a function of x 1 and x 2 , calling this function 


2 



we have 


d* 

dx 1 


u 


l 


d<a 

!h?' 


and u 2 = 


( 1 . 2 ) 


This is the approach Lagrange (l) used to introduce the stream function for two dimen- 
sional plane flows; and, later, by a similar approach Stokes (2) introduced the stream 
function for axisymmetric flows. One cam also introduce the stream function for two di- 
mensional flows by the following, slightly different, approach. A general solution of the 
continuity equation (1.1) is given by an arbitrary, as yet undetermined, function such 
that u 1 and u 2 axe related to ’F by (1.2). This approach to introduce the stream function 
for two dimensional flows can be immediatly extended to three dimensional flows. The 
extension is baised on a theorem by Jacobi (quoted in Clebsch) , which for our purposes can 
be stated as follows: The equation, 


du 1 du 2 du n 

1 1 - 1 

dx 1 dx 2 dx n 


(1.3) 


has a general solution given by (re-1) arbitrary, as yet undetermined, functions 


\f ,1 (x 1 ,x 2 ,...,x n ), 1 i f2 (x 1 ,x 2 ,...,x n ), •••, VP” 1 (x 1 ,z 2 ,...,x n ) 


(1.4) 


with u* given by its cofactor in the matrix 


/ 

u x 

U 

... U 

\ 


av 1 

av 1 

a * 1 



dx 1 

dx 2 

a*" 



a * 3 

dv 2 

a * 2 



a * 1 

dz 3 

dx n 

• • 


1 

av n ~ 1 



i 

\ 

dx 1 

dx 2 

dx n 

j 


(1.5) 


Following this line of approach Clebsch (9) introduced stream functions for three dimen- 
sional flows. We denote the two arbitrary functions for the case of three dimensional flows 
by 5 and H (instead of and ^ 2 ) and obtain from (1.5) for the velocity vector v 


v = grad'd, x grad H 


( 1 . 6 ) 


From (1.6) we note that v is normal to grad'd, thus the surface defined by a equal to a 
constant contains streamlines and a is appropriately called a stream function. We will 
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denote by £ the value of the function 5 along a given stream surface. With a similar 
discussion of H, we write our stream function equations as, 

t = Z(x,y,z)] 77 = H (x,y,z) (1.7a, 6) 

In Eqs.(l.7) the space coordinates are denoted by ( x,y,z ) to establish continuity with 
what follows latter. Another approach to introduce stream functions for three dimensional 
flows, is presented by Yih (12). Following Yih we integrate the equations of a stream line 

dx 1 _ dx 2 _ dx 3 

ny ~ = t 1 - 8 ) 

and obtain (1.7) as the integral surfaces of (1.8). Relation (1.6) is then obtained from the 
argument that since 3 and H describe stream surfaces, their gradients are normal to the 
velocity vector. The preceding discussion, for the sake of simplicity, is given for the case 
of incompressible flows. The case of compressible flows follows similarly with u l in the 
preceding discussion replaced by pu x . 

The present work takes as its starting point the existence of two stream functions that 
will describe three dimensional steady flows. From that point on we develop techniques 
for computing three dimensional flows using two stream functions. The stream surfaces in 
the present work are defined parametrically by equations such as 

* = X> y = Y[x,Z,n), and * = Z(x,£,r?), (1.9a, b, c) 

For a given value of x Eqs. (l.7a,b) and Eqs. (1.9b,c) are inverse relations of each other. In 
the present work, an important compliment to the use of two stream functions to describe 
three dimensional flows is the choice of the independent and the dependent variables used to 
describe the flow. As discussed in detail in the next section, the independent variables used 
to describe the flow are x> £> and *?; and the the dependent variables, U, the streamwise 
velocity, and Y, and Z. With these variables we studied in Ref. 13 two dimensional (plane 
and axisymmetric) parabolized viscous flows, and in Ref. 14 parabolized three dimensional 
flows through straight rectangular ducts. For such simple flows equations for U and Y, 
for two dimensional flows, and for U, Y and Z for three dimensional flows, can be easily 
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obtained by partitioning the flow into a number of appropriately defined stream tubes, and 
then applying conservation principles directly to the flow through the individual stream 
tubes. In Refs. 15 and 16 we studied two and three dimensional potential flows in term 
of these variables. In these studies equations for Y and Z were obtained by projecting the 
streamline motion on to the x-y and x-z planes. In Ref. 16 equations for Y and Z were also 
derived for the three dimensional potential flows by setting, with the help of differential 
geometry, the vorticity around a closed contour drawn on a stream surface equal to zero. 

In the present paper we present, with the help of tensor calculus, a general theory for 
studying three dimensional viscous flows using the aforementioned variables. The present 
work is restricted to steady flows. 
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II. INDEPENDENT AND DEPENDENT VARIABLES 


In this section we introduce the independent and the dependent variables. The inde- 
pendent variables are x> £, and 77. Variable x is a spatial coordinate along the main flow 
direction. Variables £ and 77 are the values of stream functions along suitably chosen two 
sets of intersecting stream surfaces. A stream surface along which the stream function 
£ is constant is referred to as a £=const. surface or as a x _r ? surface. Similar nomen- 
clature is used for stream surfaces along which 77 is constant. As discussed in the earlier 
studies (as, for example, in Maeder and Wood (11), and Yih (12)), in general, for three 
dimensional flows there is no unique way to define two sets of interacting stream surfaces. 
For a given flow there tire numerous choices for the 5 and H stream surfaces that will 
cover the given flow. However, to take advantage of the (x, £, 77) choice of the independent 
variables the general shapes of the 5 and H stream surfaces are to be selected to facilitate 
the imposition of the required boundary conditions. In the present paper we present three 
different combinations of two basic types of S and H stream surfaces that should cover 
many flows of practical interest. The two basic types of stream surfaces considered are: 
plane stream surfaces and cylinderical stream surfaces. The plane stream surfaces are not 
necessarily flat. The boundaries of plane stream surfaces intersect the flow boundaries. 
The cylinderical stream surfaces are not necessarily straight circular cylinders. The cylin- 
derical stream surfaces are nested within each other. From these two basic types of stream 
surfaces we form three different combinations, each consisting of two sets of intersecting 
stream surfaces, to model three different types of flows. These three different types of flows 
are named: (i) Plane Flows, (ii) Axial Flows, and (iii) Circulating Flows. ‘Plane’ flows are 
modeled with one set of £=constant plane stream surfaces and one set of ?7=constant plane 
stream surfaces as shown in Fig. la. This type of modeling is proposed for studying flows 
that are bounded by flat boundaries. ‘Axial’ flows are modeled with one set of £=constant 
cylinderical stream surfaces and one set of T7=constant plane stream surfaces such that 
one edge of all the 77=constant stream surfaces meet in the axis of the flow as shown in 
Fig. lb. This type of modeling is proposed for studying flows that are bounded by curved 
boundaries. ‘Circulating’ flows are modeled with one set of £=constant cylinderical 
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Axial Flows. 


8 


stream surfaces and one set of r;=constant plane stream surfaces as shown in Fig.lc. The 
intersection of the £=constant stream surfaces and r/=constant stream surfaces are closed 
curves. This type of modeling is proposed for the study of circulating flows. 

The stream surfaces are defined parametrically by the following equations: 


Plane Flows: 

x = x, y = Y{x,Z,n), 

and z 

= Z{x,t,T)), 

(2.1a) 

Axial Flows: 

x = x, r = R(x,(,v), 

and 0 

= ©(x, £>»?), 

(2.16) 

Circulating Flows: 

x = X(x, {,*}), r = R(x,£,ri), 

and 0 = x 

(2.1c) 


In Eqs.(2.1) x, y, and z are the rectangular Cartesian coordinates, and x, r, and 0 are polar 
cylinderical coordinates. For a given value of rj Eqs.(2.1) define X“£ stream surfaces with 
X and £ as the parameters. Similarly, for a given value of £ these equations parametrically 
define X'V stream surfaces with x and r? as the parameters. 

Let U(x,£,rj) represent the streamwise velocity. The dependent variables used to 
describe the flow are: 

Plane Flows: 17, Y, Z (2.2 a) 

Axial Flows: U, R, Q (2.2 b) 

Circulating Flows: U, X, R (2.2c) 

Let g x and g£ represent the coordinate vectors of the x-£ stream surfaces, and g x 
and g,, the coordinate vectors of the x-q stream surfaces (these vectors are calculated in 
the next section). Since the coordinate vector g x is common to both the x-£ and the x~V 
stream surfaces, it is tangent to the stream line defined by the intersection of the x-£ and 
the x _ *7 stream surfaces. Since U is the streamwise velocity, we note, for later use, that U 
is along g x . 
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III. METRICS OF THE STREAM WISE COORDINATE SYSTEMS 


The independent variables x, and rj form a streamwise curvilinear coordinate system. 
It is called streamwise, since, by virtue of the x coordinate lines it is aligned with the 
streamlines. In this section we present the metric of the streamwise curvilinear coordinate 
system associated with each of the three different types of flows introduced in the previous 
section. 

We start with a brief summary of the summation and the tensor notation used in this 
paper. Latin letters i, j, k, and m are used for free and dummy indices. Whenever the 
same Latin letter i, j, k, or m appears in a product, once as a subscript and once as a 
superscript, it is understood that this means a sum over all terms; thus, for example, 

3 

u% 9ii ~ 5Z = uX 9xi + + « n «7r,y. (3.1) 

«=i 

Indices x, Ct and rj represent definite directions; su mm ation is never intended on index x> 
£, or rj no matter how they appear. A subscript or a superscript proceeded by a comma 
denotes ordinary partial derivative. Thus, for example, 


r.< = 


dY 

dt 


and 


9ij,k — 


_ d 9ij 


dx k 


(3.2) 


We now return to the calculation of metrics of the streamwise curvilinear coordinate 
systems. Covariant base vectors gj of the (Xj £, r?) coordinate system are calculated from 
the transformation formula, 


Si = 


dx 3 

dx % 


S 3 


(3.3) 


where z* represent the coordinate lines of the {x,Z,v) system and x 3 are the coordinate 
lines of the (z, y, z) system for the Plane Flows, and the coordinate lines of the (z, r, 9) 
coordinate system for the Axial and the Circulating Flows. In Eq.(3.3) g y are the base 
vectors of the x 3 system and are for the (z, y, z) system 


4. = (1,0,0); 4« = (0,1,0); g, = (0,0,1) 


(3.4) 


and for the (z, r, 9) system 

g z = (1, 0, 0) ; g r = (0, cosQ, sinQ ) ; g* = (0, -RsinQ, RcosQ) (3.5) 
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Carrying out the transformation (3.3), as shown in Appendix A.l, we obtain for the co- 
variant base vectors of the Plane Flow (x, f , *?) coordinate system 

Sx = (l> ^,x > Z,x) > = (®> Kt > Z,() ; g t) = (0, Y ifl , Z )Tt ) (3.6) 

All these three covariant base vectors are expressed by the single equation, 

Si = (&x> ^.i > Z,i) (3.7) 

where the Kronecker delta, 6i x , is equal to one if » = x or zero otherwise. The covariant 
components of the metric tensor, g t y, defined by y t y = g,-.gy (where the dot between the 
vectors indicates scaler product) are given by: 

(3.8) 

(Another derivation of < 7 ,y is given in Appendix A.3) 

The determinant, g, of the matrix of the metric coefficients is 

9 = (g x -g« x g„) 2 = (y, f z„ - y,, z, ( )\ (3.9) 

where x denotes the cross product. Since yjg appears often later in the transport equations, 
we denote yfg by D, thus 

(3.10) 

To facilitate the representation of the contravariant components of the metric tensor 
we introduce Y ’*, the contravariant components of the gradient of Y, (also shown in the 
equation are Y t i, the covariant components of the gradient of Y), 

grad Y = Y'gt = Y*g x +Y’ i g ( + V’g, 

= y, S ' = Y, x g* + Y, ( gt + Y„g* (3.11) 

Gradient of Z is expressed similarly. Relations between Y •*, F,-, Z ,x and Z, are given in 
Appendix A. 2. The contravariant base vectors defined by, 


D=(y f z„-y, z, ( ). 


9ii = g.-.g,- = Mix + y.y, + z,iZj 





(3.12) 


are expressed in terms of Y ,x and Z' x and are, 

g‘ = (£ , Y>‘ , Z-‘) where 

9x = 6 x- Y >xY' i -Z,xZ' i (3.13) 

The Kronecker delta, 8 X X , is equal to one if * = x or zero otherwise. The contravariant 
components of the metric tensor, g t] \ are 

g ij = 44 + Y ’ ty,J + z,i Z^~ (3.14) 

Expressions for the contravariant base vectors and the contravariant components of the 


metric tensor in terms of Yj and Z t i are given in Appendix A.2. 

Corresponding results for the Axial and Circulating flows are: 

AXIAL FLOWS: 

gi = (6, x , R^cosQ — Q'iRsinO , RjsinO + Q^RcosQ) (3.15) 

Si] = ^ix^ix R ,i R ,j d - R 2 Q,*Q,} (3.16) 

p = R(R,j e,, - R, n e,e) | (3.17) 

g* = (g^. , R ,x cosQ — O ,x RsinO , R’ x sinQ + Q ,x RcosO) where (3.18) 

g i x = S x x ~R tX R’ i -R 2 Q a e' i (3.19) 

g 1 * = g { x g 3 x + R' { R' j + R 2 0’ l '©’ y (3.20) 

CIRCULATING FLOWS: 

g, = (X,i , R t iCosO — 6i X RsinQ , jR |t *sm© + 8i x RcosQ) (3.21) 

gij = X ti Xj + R,iR,j + 8i x 6j x R 2 (3.22) 

D = R{X, € R,„-X, n R,t) | (3.23) 

g* = ( X ’* , R’'cosO - g x Rstn 0 , R' l sinQ + g x Rcos 0) where (3-24) 

£ = (*J - X ,x X,i ~ R ,x R,i )/ R2 ( 3 - 25 ) 
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g ij = X'X'i + RJR* + R 2 g' x g J x (3.26) 

Definitions of the Christoffel symbols and Ty*, used in this study are, 

r % = g km T ijm , (3.27 a) 

2r ijk = gik,j + gjk,i - g<j,k , (3.276) 

and, in particular 

= 2?x*.x — 9xx,' i 2I\ XX = <7xx»* * (3.27c, d ) 
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IV. TRANSPORT EQUATIONS 


Equations for the transport of mass, momentum, and energy in a general curvilinear 
coordinate system have been derived previously and are given in books and review articles, 
among others, by Aris(l7), Flugge(18), and Serrin(l9). In this section we adapt these 
equations to the streamwise curvilinear coordinate system. The transport equations are 
presented in terms of the metric of a streamwise curvilinear coordinate system; and, thus, 
are equally valid for all the three curvilinear coordinate systems introduced in the last 
section. 

As discussed in section II, the velocity vector at any point is directed along g x . Thus, 
the contravariant velocity components and u v are zero everywhere. To calculate u x we 
use the relation 

U 2 = u x u 3 gij = u x u x g xx . (4.1) 

Thus, u*, the contravariant velocity components axe: 

(4.2) 

The covariant velocity components, are calculated from 

«* = 9ij = g ix u x , (4.3a) 


= «* = 0; =0. 


and are: 

(4.36) 


The contravariant velocity components, u J , in the x 3 coordinate system can now be calcu- 


u, 


_ 9xx 


9xt 


V*xx ^ XX 


U; Ur, = 


9xv 


V9. 


U . 


xx 


lated from 


u } = 


dx 3 


u* = 


dx 3 


-u* 


dx % dx 

and the physical velocity components, u(j), in the x 3 coordinate system from 


(4.4a) 


« 0 ‘) = y/iTj & 


(no sum over j) 


(4.46) 


where gjj , the covariant components of the metric tensor of the x 3 coordinate system, are 
for the (x,y,z) system 

9xx = lj 9yy ~ 1) 9zz — 1 (4.4c) 


15 





and for the (x, r, 6) system 


9xx — 1) 9rr — 1) 968 — ^ — R 


(4.4d) 


From (4.4a,b,c,d) we obtain the physical velocity components for the (see Appendix B.l): 

(4.4e) 


TJ TTY TT 7 

Plane Flows: u(x) = — — u(y) = - — ; and u(z) = -~~ 

\j9xx \/9xx \/Sxx 


Axial Flows: u(x) = - ■ ; u(r) = ; and u (6) = (4.4 /) 

\/9xx \/9xx y/9xx 

UX UR RU 

Circulating Flows: u(x) = - u(r) = ■ ■ ’* ; and u(0) = — - — (4.4g) 

y/9xx \/9xx V g xx 

The expressions for the covariant derivatives u'\j and Uj|y in terms of the Christoffel 

symbols used latter are (see Appendix B.2): 


«*'ly = u \j + u * Fyfc = + t* x T) x . 

u %\j = u*- g k i + u k T kji = u*. g xi + u x T xji 


u,' j — Ui u r tj jt — ,j u x r 




(4.5a) 

(4.56) 

(4.5c) 


In Eqs.(4.5) we have evaluated the sums over k using the fact that u x is the only nonzero 
contravariant velocity component. 

We next introduce expressions for the vorticity w, 


k = 


(Jj = e J u 


3 1» 


( 4 . 6 ) 


where the permutation tensor, is equal to +1/D , — 1 / D , or 0 depending on whether 
i,j,k is a cyclic, an anticyclic, or an acyclic sequence. From (4.6) we obtain for the 
individual components (see Appendix B.3), 




u 


i _ 


— ( U X.T Ul 0 ,x) i D 


w" = (“«.x - U x.*.)/ D 


(4.7a) 

(4.76) 

(4.7c) 
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We are now ready to develop the transport equations for a streamwise curvilinear 
coordinate system 
CONTINUITY EQUATION: 

dp 

— +div(pv) = 0, (4.8) 

where p is density. In streamwise computation of flow fields, u x is the only non-zero 
contravariant velocity component and the preceding equation becomes, 

!t + 55J (D '’“*' = 0 ' (4 - 9) 

For steady flow p t t = 0, and we obtain on integrating the remaining part of (4.9), 

Dpu x = a constant along x • (4.10) 


We select the numerical values of our stream functions such that the constant of integration 
is unity. The steady state continuity equation thus becomes, 



(4.11a) 


and with the help of (4.2), 

U = * (4.116) 

Relationship of Eqs.(4.11b) and (4.4e) of this section with Eq.(1.6) of the introduction 
section is discussed in Appendix C. 


MOMENTUM EQUATION: 


D\ dv ,• _ 

"Si = '’ar + ' >uv ..' = F ’ 


(4.12) 


where F is the force (surface and body) per unit volume. Once again, since u x is the only 
non-zero contravariant velocity component, the covariant and the contravariant momentum 
equations become, respectively, 


P^i ( U,flr *) + PUX = F ' 9 ' ' (4.13a) 

P§- t {» X 9x) + P» X \x9i = . (4.136) 
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For steady flow Eqs.(4.13) become, 


pu x m\ x = Fi 

(4.14a) 

pu x u*| x = F x 

(4.146) 


Development of the covariant components of the momentum equation: 
With the help of (4.5b) we rewrite Eq.(4.14a) as 

pu x (u x x g x i + u x r xx ,) = Fi . 

With the help of continuity equation (4.11a), (4.15) becomes 

1 d_ f 1 \ 1 

Dd X \ P D) 9xi+ pD* Txxi ' Fi ' 

Expressing via (3.27c) we obtain after some algebra, 


(4.15) 


( 4 . 16 ) 


pDg x i iX g xi ( pD ) tX ^ P^9xx — P* D 3 F 


(4.17) 


Another variation of momentum equation (4-14a) is obtained by using the expression for 
u,-|y given in (4.5c), 

pu x {m lX — u x I\ xx ) = F x . (4-18) 

Using (3.27d) to express I\ xx and that 

u x>* = (y/9xx^),i ~ V^xx^.i + -jUg XXj i/ y/g xx = V?xx^.» + T^* 5 **’* 


we obtain from (4.18), 


pu x {u itX - u Xti + Vfj xx U,i) = Fi . 


(4.20) 


With the help of the continuity equation and some rearrangement we obtain from (4.20) 


pUU t i — F{ £}( u »,x u x>») • 


(4.21a) 


Writing out the individual components, and using the expressions for u given in (4.7), we 
have, 

(4.216, c,d) 


pUU tX = F x , pUU t t = F{ - w", pUU <t7 = F n + 
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Development of the contravariant components of the momentum equation: 

With the help of (4.5a) we rewrite Eq.(4.14b) as 

pu*(u‘ x + u*ry = F*. (4.22) 

With the help of continuity equation (4.22) becomes 

i-(JL)=.DF*, I i x = ?D*Ft, rj x = pD'F* . (4.23a, b,c) 

ENERGY EQUATION: 

— + div(pev) = ~div q + dt't;(T.v) + pv.b , (4.24) 

£ 

where e is the energy (sum of internal, kinetic, and potential) per unit mass, q the heat 
flux, T the stress tensor, and b the body force per unit mass. Once again using the fact 
that u x is the only non-zero velocity component, we obtain from (4.24) for steady flow, 

i A( W) - . | (4.25) 

CALCULATION OF F { : 

We separate the pressure, the viscous, and the body force contributions to Fi and write 


Ft = —p,i + fi + phi 


(4.26) 


where 6,- is the body force, and the viscous contribution, 

fk = Ski r = Ski r% = Ski [i A + r> m r;. m ] (4.27) 

with the viscous stress tensor r t} given by 

r‘> = -l^Us* + n(s im u‘L + S im u *\ m ) • (4-28) 

The stress tensor T' : is given by 

r'i = -pg ij + T ij (4.29) 
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V. POTENTIAL FLOW 


In this section we discuss the inviscid approximation to the transport equations pre- 
sented in the preceding section. The main purpose of this discussion is to illustrate the 
procedure for obtaining transport equations for a particular type of flow (for example, 
Plane, Axial, or Circulating) from the general transport equations. 

For inviscid flow and neglecting body forces, 


Fi = ~P,i , 


(5.1) 


and the momentum equation (4.21b) becomes, 

P uv ,x = -p, x . 


(5.2) 


We further assume that the flow is isentropic, and thus (see Appendix D.l), 

where 7 is the ratio of specific heats. With the help of (5.3) we integrate (5.2) and obtain, 

(5.4) 


U 2 IP 

— — I = const. 

2 q-lp 


With the help of the continuity equation (4.11b), we rewrite (5.4) as, 


9xx .'IP A 

1 = const. 


2p 2 D 2 7 — 1 p 


(5.5) 


From (5.3) and (5.5) we obtain, 


P,i — P{9xx/lp 2 D 2 ),i 


(5.6) 


From here on subscript i stands only for £ and rj. From Eqs.(4.17), (5.1), and (5.6) we 
obtain, 

pDg x i <x ~ 9 x'(pD),x ~ 2 pDQxx,* = P 3 ^ 3 {9xxl‘^P ^ ),t (5-7) 

With a little algebra (5.7) can be rewritten as (see Appendix D.2), 


{9x'/pD),x {.9xxlpF>),i — 0. 


(5.8) 
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Equation (5.8) can also be obtained by setting the vorticity w* (see Eqs. 4.7b and 4.7c) 
equal to zero and using the continuity equation (4.11a) along with Eq (4.3a) to express u,-, 
the covariant components of the velocity. 

To separate the effects of compressibility we rewrite (5.7) as, 

9x*,xD ~ 9x*D,x ~ 9 XX ,%D + g xx D t i = — ( 9x*P,x ~ 9 XX P , «') ( 5 - 9 ) 

For incompressible flows the right hand side of equation (5.9) vanishes and we have, 

9x*,xD — g X {D tX — g xx ,\D + g xx Di = 0 (5.10) 

Equation (5.8) is our basic equation to describe isentropic flows. Upon substitution of 
appropriate expressions for g x i, g xx , and D from section III for a given type of flow, Eq. 
(5.8) will yield the streamwise isentropic flow equations for that particular flow. In the 
following we illustrate the procedure for obtaining isentropic flow equations for a particular 
type of flow from the general equation (5.8). 

Incompressible Plane Flows: Upon substituting derivatives of <?,y and D from Eqs. (3.9) 
and (3.10) into (5.10) we get the following equation for the x-£ stream surfaces, 

D{Y,xxY,Z + Z,xxZ,z ~ Y tX ^Y jX — Z jX ^Z tX ) — g x ^(Y iX ^Z >r) + Y^Z tXri — Y iXV Z t £ + Y jV Z tX ^) 

+9xxiX,zzZ,v + Y t Atn ~ W.« + = °> ( 5 - lla ) 

and the following equation for the x _r ? stream surfaces, 

D(Y xx Y >r , + Z, xx Z tTI — Y tXn Y tX — Z, xv Z tX ) — g xrt (Y tX ^Z >ri + Y t ^Z iXn — Y iXn Z } { — Y )V Z tX $) 

+ 9x xiy.tnZ,n + Y ti Z, m - Y„ n Z ti - Y.Z^) = 0, (5.116) 

where D, g x $, and g xr] are given by Eqs. (3. 9) and (3.10). We have two equations for two 
unknowns Y and Z. These equations were solved numerically in Ref. 16 for (i) flow through 
a rectangular diffuser with an offset and a change in the aspect ratio, and (ii) flow through 
a duct whose cross-section changes from a square to a rhombus. For a two-dimensional 
flow in the plane 

Y t n = Z, x = Z,t = 0, and Z tV = 1 
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and Eq.(5.11a) reduces to, 

KtKxx - + (1 + r’)y,« = o (5.12) 

This equation was solved for flow past an infinite cylinder placed symmetrically between 
two parallel plates in Ref. 15. 

Equations for axial and circulating flows are obtained in a similar fashion. Solution 
of the axial flow equations for flow through superelliptic transitional ducts is given in a 
forthcoming article. We next present one special case each of axial and circulating flow. 
AXIAL FLOW SPECIAL CASE: For axisymmetric flow in X <, 

= 0 iX = R iT) = 0, and 0 )T) = 1 (5.13a) 

and we have from Eqs.(3.14) and (3.15) 

9xx = 1 + R,x* 9xt = ***.«, an ^ ^ = RR>( (5.136) 

With the help of Eqs.(5.13), Eq.(5.9) becomes, 

R,(R,xx ~ ^RixR(R>x^ d" d" R,t/R = {R,t/ p)IR'xR>£PiX ~ (5.14) 

Equation (5.14) can be used to study isentropic flow through nozzles. 

CIRCULATING FLOW SPECIAL CASE: For two-dimensional circulating flows in 
x-R 

X, x = X t £ = R t r, = 0, and X tfl = 1 (5.15a) 

and we have from Eqs.(3.22) and (3.23) 

g xx = R 2 + R 2 x , g X £ = R, x R,z, and D = —RR,z (5.156) 

With the help of Eqs.(5.15), Eq.(5.10) becomes, 

R%R,xx ~ 2 R,xR,tR,xt + (R + R,x)R,a ~ RR,z = ® ( 5 . 16 ) 

We note that the free line vortex given by 

R = const, e* (5-17) 
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is a solution of Eq.(5.16). To see the relation between (5.17) and the more familiar form 
of the line vortex equation, we obtain from Eqs.(4.4g) for the solution given in (5.17) 

u(x) = 0, u(r) = 0, and u(6) = U (5.18) 

For compressible flows £ and rj are related to the mass flow rate with the SI units of \/kg/s. 
For flows with constant density we rewrite pD for the circulating flow as 

pD = R(X t U /jp) R( n /^p) - R,{i/^p)) (5.19) 

For such flows we absorb y/p into the definitins of £ and rj, which now become related to 
the volume flow rate with the SI units of >/m 3 /s. Thus, for flows with constant density 
the continuity equation (4.11) becomes in terms of the newly defined £ and rj, 

U = (5.20) 

By the use of (5-15b) and (5-17) we obtain from (5.20) 

77 _ _ _ R _ 

D RR,z ~ 

From (5.18) and (5.21) we obtain the familiar form 


const. R 


(5.21) 


of the line vortex equation. 
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VI. PARABOLIZED FLOWS 


In this section we present an approximation to the viscous forces which for two dimen- 
sional flows reduces to the well known boundary layer approximation. For lack of a better 
name we have called it parabolized flow approximation. 

We assume that all contributions to the viscous forces that arise from the field variation 
of the curvature of the coordinate system (that is, from the terms involving Christoffel 
symbols) are negligible. Thus, (a) 

and, (b) u*|y are equal to u*y which have only one non zero component a*. Approximation 
(b) leads to that 


r xx = , T x( = T Sx = , 

r x ” = r" x = f ,g" m u x m , r« = (-2/3) W «u* , t" = (-2/3)« S ’"'«* , (6.2) 

and that is negligible. We also assume that the contributions of the cross-stream 
velocity derivative terms and u* to the stress dominate over the contribution of the 
u* terms. Thus, the only non-negligible terms of the stress tensor we axe left with are, 


T ™ = »(gxtu^ + g xr, ux) 
r x5 = T Zx — ^(g££ u * + g^u*,) 

T XV = T VX = M (g^ U X + g^u*) 


(6.3 a) 

(6.36) 

(6.3c) 


We further assume that in Eq.(6.1) the contribution of the term involving the derivative 
with respect to x is negligible, thus 


P = — — (Dr 1 ' 5 ) + — — (Dr*') 
1 Ddt X } D dr) V ' 


(6.4) 


From (6.3) and (6.4) we have that 

fX = hh V D + •‘"•s)] + V D (•"‘“S +ff ” u '0] < 6 - 5 > 


24 



and that ft and f v are negligible. For the covariant component of the viscous force we 
now have 

ft = Qijfi — 9i x f X » (6.6a) 

which leads to 

fx = 9 X xf X » A = 0£x/ x > fn = Srj*/ X • ( 6 - 66 > c » rf ) 


We conclude this section by deriving the viscous force term used in the boundary layer 
approximation for two dimensional flows from the viscous force approximation given by 
Eq. (6.5). From equation (6.5) we obtain for two dimensional plane flows in (x,y) with 
X = 2 , £ = y, g vy = 1 , D = 1 , and the physical velocity along x, u = y/g xx «* = « z , 



(6.7a) 


and for axisymmetric flows in (x, r) with x = £ = r, g rr = 1, D = r, and the physical 

velocity along x, u = y/g^u x = u x , 



(6.76) 
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VII. CONCLUDING REMARKS 


A theory to compute three-dimensional flows using two stream functions is presented. 
The values of two stream functions along with a spatial coordinate are used as independent 
variables. Since the value of a stream function is constant along the solid boundaries, 
this choice of variables makes it easy to satisfy the boundary conditions. The dependent 
variables employed axe the streamwise velocity and two functions that define the stream 
surfaces. Thus, the method generates a boundary fitted grid that is aligned with the flow 
streamlines as part of its solution. 

For three-dimensional flows so far there is no general theory to define two sets of 
intersecting stream surfaces to cover a given flow. In the present work we have presented 
three different combinations, each consisting of two sets of intersecting stream surfaces, 
which should cover many flows of practical interest. However, a general theory regarding 
the selection of two sets of intersecting stream surfaces that will cover a given flow for the 
most efficient computation of the flow is needed. 
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APPENDIX A 


A.l In this appendix we derive the covariant base vectors of the streamwise curvilinear 
coordinate system, x*, for the three different types of flows introduced in section II. 
PLANE FLOWS: 

Base vectors, g,-, of the cartesian coordinate system x* = (x, y, z) 

g x = (1,0,0); g y = (0,1,0); g z = (0,0,1) (A.1.1) 


are transformed by 

dx^ 

Si = jjrfc (A.l. 2) 

to obtain the base vectors, g t -, of the x* = (x. £,w) coordinate system. From Eqs. (A. 1.1) 
and (A. 1.2) we obtain, 


dx . . dy , , dz . . . dx 

g. = gjd.o.o) + ^(0,1,0) + ^(0.0,1) = 


dy dz 
dx 1 ' dx 1 


(A. 1.3) 


Now using the transformation equations 


* = X, y = P(x,f.>7) } and z = Z{x,^v), (A.1.4) 

we obtain Eqs. (3.6) from Eqs. (A. 1.3). 

AXIAL FLOWS: 

Base vectors, g,-, of the polar cylindrical coordinate system x* = (x,r, 6) 

g x = (1,0,0); g r = (0 ,cos0,sind) gg = (0,~r sind,r cosd) (A.1.5) 

are transformed by Eq.(A.1.2) to obtain the base vectors, g», of the x* = (x, £,r?) coordinate 
system, 


dx , .dr , . . /.x d0 . . 

g, = + ^7(°’ cos0 > 5m0 ) + ^7(O,-rstn0,rcos0) 

.dx dr dd . dr . dd . 

= (-r— • , -r—^cosd - —~^rsin0 , ~-stw - — T r COSO). 


' dx x dx x dx x 

Now using the transformation equations 


(A. 1.6) 


x = x, r = R{x,Z,ri), and 0 = Q(x,Z,v), 


(A. 1.7) 
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we obtain Eqs.(3.15) from Eqs.(A.1.6). 

CIRCULATING FLOWS: 

Base vectors, g,, of the polar cylindrical coordinate system z* = (z, r, 6) given in (A. 1.5) 
are transformed by Eq.(A.1.2) to obtain the base vectors, g,-, of the z* = (x, £, r?) coordinate 
system, once again, given by Eq.(A.1.6). Now using the transformation equations 

* = ^(X»£,*7), *■ = R{X,Z,V), and 6 = x (A. 1.8) 

we obtain Eqs.(3.21) from Eqs.(A.1.6). 


A. 2 In this appendix we present the contravariant base vectors and the contravariant 
components of the the metric tensor in terms F t - and Z ,. 

PLANE FLOWS: 


g x = iiill = (1,0,0) 

£ — * gx _ W,r\Z,X ~ ^.X^V? » Z,t) » ~Y,rt) 

8 D D 

n = Sx x g$ = ( V .x z .( ~ Y ,i z ,x » ~ z ,i i 
8 D D 


(A. 2. la, 6,c) 


The contravariant components of the metric tensor, g x * , defined by g tJ = g‘.g J , are: 


jXX = l 

(A. 2.2a) 

S x( =S ix = [Y, v Z, x -Y, x Z,„)!D 

(A. 2.26) 

9 X ” = 9 nx = (V,x Z,i - y,( Z,x)/D 

(A. 2.2c) 

9 (( = („*£)’ + (y» + Z’ )/.D 2 

(A.2.2d) 

S S " = „■>« = gXlgXV _ (y £ y„ + Z{ z,„)/D» 

(A.2.2e) 

S’" 1 = (s*’) 2 + (rf + Z? ( )/D 2 

(A. 2.2/) 


Relations between the covariant and the contravariant components of gradY and gradZ 
are obtained from 

Y’* = g ij Y ih Z ,x = g i] Zj, (A.2.2g) 
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with g'J given in (A2.2), and are, 

Y ' x ~ 0 , Y>t = Z ifl /D, Y" = -Z ti /D 

Z ,x = 0 , Z’t = -YJD, Z* = Yt/D 
Y’% = Z'% = 1 and y*% = Z'% = 0 

AXIAL FLOWS: 

S x = (1,0,0) 

~ &l R ( R ''> e .x - R ,x Q ,r>) 7 {Rr,sinO + e >ri Rcos0 ) , {-R tri cos0 + e t r,RsinQ ) j 

S ~ £ W^.x 0 ^ ~ R ,t e ,x) , ~(R^sin$ + 9 ti RcosO ) , (R^cosd - O^RsinQ)] 

(A.2.3o, b, c) 


gXX = 1 ( A . 2 . 4 a ) 

j*< = g ( * = 0,* - a, x e,,)/p {A 2 ib) 

a *' = = R ( R, X ©, ( - *, £ 0 ,*)/b (a. 2 . 4 C ) 

9 £{ = ( 9 «) 2 + ( J ? 2 + ft 2 © 2 ,)//) 2 ( A . 2 . 4 d ) 

!,£■> = *■>< = S* V - (ie >£ B , + i? 2 0, f e,,)/z> 2 (A.2.4e) 

9”’ = (0 2 + (fi, 2 { +fi 2 ef £ )/Z) 2 (A. 2 . 4 /) 

*'* = °. R' ( = R&,JD R'" = -RG t /D ( A .2Ag) 

0' X = °. Re> l( = -R,,/D, Re-"=R A /D (A.2Ah) 

R,<R ,< = je 2 ©'-0, ( = l an d fl-‘jJ0 j£ = RQ.‘ Ri = 0 (A . 2 . 4>) 


(A.2.2h) 

(A.2.2i) 

(A.2.2j) 
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CIRCULATING FLOWS: 


g x = (0, — sinQ,cosQ)/R 

> (-^,» 7 -^iX — RMsinQ — X ifl RcosQ, [X tf) R )X — X tX R tTI )cosQ — i2.A jrf stn©] 

g* 7 = {R,xX,t ~ R ti X a )sinQ + X^RcosQ, (X tX R t ^ — X^R iX )cosQ + X^RsinQ] 

(A. 2.5a, b, c) 

g xx = 1/R 2 (A. 2.6a) 

gXt = fix = (X, B 1X - X, x *,„)/*!> (A.2.64) 

9*" = s’* = (X, x R, ( - X,< JW/M (A. 2.6c) 


9« = B 2 [(s* £ ) 2 + (X 2 + Jif„)/X> 2 ] (A.2.6d) 

g in = ft = B 2 [(sXC s « _ (X C X, + * >{ JE.,)/Z? 2 ] (A.2.6e) 

S ’” = ^ 2 ((?*") 2 + ( X 2 { + X 2 f )/ n 2 ] ( A . 2 . 6 /) 


X-* = 0, X’ $ = RR,r,/D, X ’•» = —RR't/D {A.2.6g) 

R ,x = 0, R-* = —RX'fj/D, R ,TI = RX't/D (A.2.6h) 

X^Xs = R''R t i = 1 and AT-’R < = R>‘X, = 0 (A.2.6t) 


A. 3 In this appendix we calculate gtj from the transformation formula, 

dx m dx n . 

9ij ~ dx •' <9x> 5mn 


(A. 3.1) 


where g mn are the covaxiant coefficients of the metric tensor of the (x,y,z) coordinate 
system for the Plane flows, and of the (x, r, 8 ) cordinate system for the Axial and the 
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Circulating flows. For all three flows (Plane, Axial, and Circulating), g mn = 0 if m ^ n 
and (A. 3.1) reduces to, 


dx 1 dx 1 . dx 2 dx 2 „ 3x 3 dx 3 

dxi 911 + ax' axi 922 + a x < ax> 933 

(A.3.2) 

PLANE FLOWS: 

By the use of following transformation equations and gu 

x = x, y = Y(x,t,v), z = Z(x,Z>ri), 

(2.1a) 

9xx = 1) Qyy = lj Qzz = 1 

(4.4c) 

we obtain from (A.3.2), 

9ij = + Y ,i Y j + Y ,iZ,j 

(3.8) 

AXIAL FLOWS: 

By the use of following transformation equations and gu 

x = x, r = R{x,£,ri), 0 = ©(x, £,r?), 

(2.16) 

9xx = 1; 9rr = 1; 9ee — r 2 = R 2 

( 4Ad) 

we obtain from (A.3.2), 

9ij = "h "h R ®,*® |J 

(3.16) 

CIRCULATING FLOWS: 

By the use of following transformation equations and gu 

X = x{x, Z,v), r = R(x,Z,ri), 8 = X 

(2.1c) 

9x x = 1; 9rr = 1; 9ee = r 2 — R 2 

(4.4<f) 

we obtain from (A.3.2), 

9ij = X,iX,j + R,i R,i + 8i x 6 jx R 2 

(3.22) 
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APPENDIX B 


B.l Equations (4.4e)-(4.4g) are obtained from Eqs. (2.1), (4.4a), (4.4b) and (4.4c). 
Recall 


Hj) = VIjj = Vljj ^ « x 


(no sum over j) (4.4a, 6) 


and 


Six 1) 9yy — lj 9zz — 1 
Six = 1; 9tt = 1; 9ee = r 2 = R 2 


(4.4c) 

(4.4d) 


The details are: 

PLANE FLOWS: 


AXIAL FLOWS: 


u(z) = y/§T x ^u x 

Hy) = 

H z ) = 


U 


y/9xx 

Y, X U 

y/9xx 

Z, X U 

y/9xx 


\ /x — x U 

U(x) = = — 

U(r) = yfct. u* = 

d X y/9 K 

pr~ de x RQ ,x U 


(B.l. la, 6, c) 


(B.1.2a, 6,c) 


CIRCULATING FLOWS: 


u(x) = y/OZ 

u(r) = 

H 9 ) = Vse* 


X,xU 

y/9 xx 

ML 

y/9xx 

RU 


y/9 xx 


(B.l. 3a, 6, c) 
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B.2 In this appendix we derive the expressions for the covariant derivatives given in 
section IV. We start with the expression for tx*jy , the i-contravariant component of the 
j- derivative of v (see Ref. 18 Eq.(S.ll)), 

uii = u‘j- + tt*r}* (£.2.1) 

This is Eq.(4.5a). On multiplying (B.2.1) with and using 

u \j9im ~ 9im)\j = u m\j j 
u j9im = Uj9km and 

^jk9int T'jkm ~ r* kjm (B.2.2fl, 6, c) 

we obtain 

«* m |y = 9km + u k T k jm {B.2. 3) 

Equation (B.2.3) is Eq.(4.5b) with the free index t replaced by m. To obtain Eq.(4.5c) we 
start with Eq.(5.13) of Ref.18, 

u t |y = u, t j Uk r ,y (B .2.3) 

On multiplying the last term on the right hand side of (B.2.4) with gkig kl = 1, and then 
using g kl to raise the index on u* and using gki to lower the index on T k } - we obtain 
Eq.(4.5c). 

B.3 From 

u> k = t iik Uj\i (4.6) 

we obtain for w", 

u" = = («elx - «xl ()/» (B.3. la) 

Similarily we obtain, 

= ( w xU ~ U M/ D [B.Z.lb) 

w x = [ u n\t ~ u z\ti)/D (B.3.1c) 

Now, 

Ui\j - Uj\i = Uj,y - tt*r f y* - + u k Tjik = Ui,j - Uj t i [B. 3.2) 

where we have used r t y* = Ty t fc. Equations (B.3.1) and (B.3. 2) lead to Eqs. (4.7) 
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APPENDIX C 


From Eqs.(4.11b) and (4.4e) we obtain 

P u i x ) = P*(v) = />«(*) = (C.l) 

For a fixed value of £, say £o> equations 

z = x, y = ^(x.£o,*?)i * = ^(x. £o,n), (C-2) 

describe the stream surface E. The gradient of E is given by the cross product of the base 
vectors of E, g x and g n , thus, 

grad a = (l ,Y >X ,Z X ) x (1 ,Y tT ,,Z v ) (C.3) 


Similarly we have, 

grad H = (1 ,Y, 0 Z ( ) x (1 ,Y, X ,Z X ) (C.4) 

From(C.l),(C.3), and(C.4) we obtain, 

pv = grade, x grad H (C.5) 

Equation(C.5) is the compressible version of (1.6). 
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APPENDIX D 


D.l We start with the well known thermodynamic relation 

Tds = dh — dp/p (D.l.l) 

where s is entropy and h is enthalpy. For isentropic flow of an ideal gas, 

ds = 0, and dh = d(c v T) — d — d ^ • (D.l. 2) 

Upon substituting relations from (D.1.2) into (D.1.1) we obtain (5.3). 

D.2 Upon dividing (5.7) by p 2 D 2 we obtain, respectively, for the first, second, and 
third term on the L.H.S. and for the term on the R.H.S., 


9x* ,x 
pD ’ 

[D. 2.1) 

p 2 D 2 (f> D ).X-9xt ( pD ) x > 

(£>.2.2) 

^ 9 xx ,* 

2 pD ’ 

)£>.2.3) 

2 (Wp d ),» - 2pD 9xx ’ x +9xx ( p d) { ' 

(D.2.4) 


Upon combining (D.2.1) with (D.2. 2) we obtain the first term on L.H.S. of (5.8), and upon 
combining (D.2.3) and (D.2.4) we obtain the second term on the L.H.S. of (5.8). 
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